The curvature perturbation in a box 
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The stochastic properties of cosmological perturbations are best defined through the Fourier 
expansion in a finite box. I discuss the reasons for tliat with reference to the curvature perturbation, 
and explore some issues arising from it. 
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I. INTRODUCTION 

Strong observational constraints on the primordial cur- 
vature perturbation Q make it an important discriminator 
between models of the very early Universe. After smooth- 
ing relevant quantities on the shortest scale of interest, 
C may be defined as the fractional perturbation 5a /a of 
the locally-defined scale factor a(x, t). Equivalcntly, it 
is the perturbation in the number N{x,t) of e-folds of 
expansion, starting on an initial 'flat' slice of spacetime 
(where Sa — 0) and ending on a slice of uniform energy 
density at time t. The spacetime threads of constant x 
are taken to be comoving. 

Any choice of the initial 'flat' slice will do, as long as 
the smoothing scale is outside the horizon at that stage, 
because the expansion going from one 'flat' slice to an- 
other is uniform. The final slice is to be located before the 
smoothing scale re-enters the horizon, and on cosmologi- 
cal scales it should be late enough that ( has settled down 
to the final time-independent value that is constrained by 
observation. With the smoothing scale outside the hori- 
zon, the evolution at each point is expected to be that 
of some uriperturbed universe (the separate universe as- 
sumption [ii,[3,ii,ii. 

During inflation, one or more of the scalar field per- 
turbations is supposed to be practically massless during 
inflation (m^ ^ H^). For these 'light' fields, the vac- 
uum fluctuation is promoted to a classical perturbation 
as each scale leaves the horizon. To calculate SN one 
smooths the light field perturbations on a comoving scale 
shorter than any of interest. A few Hubble times after 
the scale leaves the horizon the light fields have classical 
perturbations. Their values at this epoch are supposed 
to determine iV(x, t) and hence Taking for simplicity 
just one field a we have 



where 



and 



C(x,i) = SNiai^)) 

= N{a{x))^N{a) 



Scr{x) = cr(x) — cr. 



dN 
da 



N" 
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and so on. 

The final expression for C is a power series in the field 
perturbation defined on the initial slice. The unper- 
turbed value of a is defined as its spatial average of each 
field, which simplifies the analysis. 

We are interested in the Fourier components of C: 



Ck 



(6) 



and similarly for 6a. The perturbation Sa^ is supposed 
to originate from the vacuum fiuctuation, and can be 
considered 0] as classical starting a few Hubble times 
after the epoch of horizon exit k = aH = a, where a 
is the scale factor. At this stage the correlators of Sau 
on scales not too far outside the horizon can easily be 
calculated using perturbative quantum field theory. As- 
suming Lorentz invariance and a quadratic kinetic term, 
da is almost gaussian. The subsequent separate-universe 
evolution of 5a {x.) then gives the correlators of 6a at the 
initial epoch. Finally, Eq. ^ gives the correlators of C 
which can be compared with observation. 

To avoid assumptions about the unknownable Universe 
very far beyond the present horizon Hq^ , this calculation 
should be done within a comoving box, whose present size 
L is not too much bigger than H^^. This situation has 
not been discussed much in the literature for two reasons. 
First, with the box size not too much bigger than Hq^ 
(minimal box) the linear term of Eq. ([3]) dominates. Then 
correlators are practically independent of the box size, so 
that only a need be specified. Second, if a is the inflaton 
in single-field slow-roll inflation and we use a minimal 
box, then a can be calculated from the inflation model. 

On the other hand, a may have nothing to do with 
inflation, as in the curvaton model 0]. Also, we should 
know how to handle the dependence on box size as a mat- 
ter of principle. In this note I look at some of the issues 
raised when one takes seriously the box size, developing 
some earlier work [Bl, 1^ . 



II. CORRELATORS OF ( 

A model of the early Universe will predict, not C(x) 
itself but correlators, {({x)({y)) and so on. The () is 
an ensemble average, which in the inflationary cosmology 
becomes a Heisenberg-picture vacuum expectation value. 
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As the vacuum is invariant under translations, so are the 
correlators. 

The vacuum is also invariant under rotations, and so 
are the correlators.^ Invariance under translations and 
rotations constrains the form of the correlators. One 
usually considers only the two, three and four-point cor- 
relators. 

The two-point correlator is 

(CkCk') = (2vr)=^5(3)(k + k')Pc(fc), (7) 

where is the spectrum. It is useful to define = 
{k^ /2tt'^)P(^, also called the spectrum. The mean-square 
of ( is 

(C^W) = £7n(fc)f , (8) 

where /cmax is the scale leaving the horizon at the initial 
epoch, and the infrared cutoff is provided by the box size. 
On cosmological scales observation gives Vi^ — {5x 10~^)^ 
with little scale dependence. 
The three-point correlator is 

(CkiCk.Ck,,) - {27Tf6^^\-k,+k2+-k3)Bc{h,k2,h), (9) 

where B,^ is the bispectrum. The connected contribution 
to the four-point correlator is 

(CkXk.CkaCkJc = (27r)35(3)(j^^ + k2 + k3 + k4)Tc. (10) 

The trispectrum is a function of six scalars, defining 
the quadrilateral formed by {ki, k2, ka, k4}. There is also 
a disconnected contribution: 

(CkiCk2Ck3Ck4)d (CkiCk2)(Ck3Ck4) + perms. (11) 

For any correlator, the connected correlator is the one 
that comes with an overall delta function. If the two- 
point correlator is the only connected one, is said to 
be Gaussian. Data are at present consistent with the 
hypothesis that ( is perfectly gaussian, but they might 
not be in the future. 

As a result of translation invariance, the ensemble aver- 
ages (((x)((y)) etc. can be regarded as a spatial average 
with fixed x— y. This is the ergodic theorem, whose proof 
for correlators is very simple. Considering the two-point 
correlator the definition of Pq gives 

(C(y)C(x + y)) = J P{k)e^--d^k. (12) 



In considering translations and rotations one takes x to denote 
Cartesian coordinates defined in the unperturbed universe. In 
the perturbed universe these same coordinates cover the curved 
spacetime slice of fixed t, in a way that is defined by the threading 
which is here taken to be comoving. 



On the other hand, the spatial average is 

L-'J C(y)C(x + y)d^y 

= L-3(2^)-6 I (;^(;^,e^[^ y+'^'<^+y)U^yd''kd'k' 

= L-^{2t:)-^ J 6^{k + k')CkCk'e'^'""d^kd^k'. (13) 

In the final expression, CkCk' can be replaced by its en- 
semble average (CkCk'}i because each volume element 
d^kd^k' can contain an arbitrary large number of of 
points. (Remember that we are working in a finite box 
so that the possible momenta form a cubic lattice in k- 
space. Within a cell, the Fourier coefficients are uncorre- 
lated because of the delta functions in Eqs. (O, ([5]), and 
(fTO|) and so on.) Writing (CkCk') in terms of P{k) and 
using the rule [j^(k - k')]^ = L^S^{'k - k'), we see that 
the spatial average is indeed equal to the ensemble 
average ([T2|) . reproduces the ensemble average, and the 
same argument works for higher correlators too. 

From this proof it is clear that the ergodic theorem 
works in a finite box, with the usual proviso that the 
box size is bigger than scales of interest. Its proof re- 
lies just on translation invariance, which makes mathe- 
matical sense in the finite box because of the periodic 
boundary condition. 

III. THE CORRELATORS OF Sa 

The correlators of 6a have the same form as those of 
C. In particular 

(^k^k') = i27TfS^^\k + k')P,{k) (14) 
(0-ki0-k2O-k3O-k4)d = (o-kiCrk2)(crk3 0-k4) + perms. (15) 

The chosen box should be well inside the horizon at the 
beginning of inflation. After a few Hubble times, the uni- 
verse inside the box is then expected [3] to become prac- 
tically homogeneous and isotropic at the classical level. 
The conversion of the vacuum fluctuation into a classical 
perturbation Sa begins when the box leaves the horizon 
at the epoch k = aH. 

We observe scales k > Hq where Hq is the present 
value of the Hubble parameter. Taking inflation to be 
almost exponential there are at most 60 or so Hubble 
times between horizon exit for the scale k — Hq and the 
end of inflation. Of most interest is the perturbation on 
cosmological scales, which leave the horizon during the 
first 10 or so Hubble times. Smaller scale perturbations 
could also be of interest, for example to form black holes 
at the end of infiation. 

We can distinguish between two kinds of box. A min- 
imal box, for which Iii{LHq) is not too big, will leave the 
horizon just a few e-folds before the observable Universe 
leaves the horizon, and hence not too many Hubble times 
before all observable scales leave the horizon. A super- 
large box on the other hand, with very large In(LiJo), 
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would leave the horizon very many e-folds before the ob- 
servable universe. We will see how to calculate things 
in a minimal box, and poin to the difficulties that are 
encountered if one considers instead a super-large box. 

With the minimal box, Sa remains small at least while 
cosmological scales leave the horizon. On the usual as- 
sumption that (T is canonically normalized, the pertur- 
bation 6a generated from the vacuum is then almost 
gaussian on scales not too far outside the horizon. Its 
spectrum is [l^l 



(16) 



where Hk is the Hubble parameter at horizon exit. The 
bispectrum 11 1 and trispectrum jT2i] of Sa have been cal- 
culated, and are suppressed by slow-roll factors. (These 
are \H/H'^\ ^ 1, needed for inflation, and parameters 
involving derivatives of the potential V{a) which have to 
be small to justify the perturbative quantum fleld theory 
calculation.) 

As we saw earlier, the correlation functions deflned by 
the spectrum, bispectrum and so on may be defined as 
spatial averages within the box. With a minimal box it 
is reasonable to assume that these are the same as the 
spatial averages within the region of size H^^ around us, 
that can actually be observed. 

If we consider instead a super-large box several difficul- 
ties arise. We have to assume that enough inflation took 
place for the box to exist, and we need to understand the 
field theory starting from the era when the box leaves the 
horizon. When observable scales leave the horizon, the 
perturbation that has already been generated on larger 
scales my be large, which would complicate the calcu- 
lation of 5a. And even when the calculation has been 
performed, the spatial averages represented by the corre- 
lators may have nothing to do with the spatial averages 
that are actually observed. 



IV. THE MEAN VALUE a 

A. With £7 the inflaton 

Now we come to the mean value ct, of the field within 
the box. If a is the inflaton in a single-field slow-roll 
infiation model, and we use a minimal box, then ct when 
N e-folds of inflation remain is given by 



N{a) = Ml r ^da. 



(17) 



This follows from dN = Hdt and the slow-roll approxi- 
mations 



Bffcr = 

2 _ 



H/H 



-V 
V 

~a\ 



(18) 
(19) 
(20) 



The inflation model will give the field (Tend at the end of 
infiation, and the post-inflationary cosmology determines 
N when a given scale leaves the horizon, with iV ~ 50 or 
so for cosmological scales in the usual cosmology. 

In this very special case, C becomes time-independent 
soon after horizon exit. As a result the small pertur- 
bation C = 5N ^ 10~^ is then given in terms of the 
potential and its derivatives [ill , fisj : 



1 l.o 



1 V 



P) -^—5a 



(21) 



where 2e = Ml(V'/Vf, rj = M^V" /V and = 
M^V"'V' /V"^ . With a minimal box the linear term dom- 
inates, and taking the initial slice to be soon after horizon 
exit on a given scale we then arrive at the famous result 



1 V H 
'mIv'2^ 



(22) 



where the right hand side is evaluated at horizon exit. 
The bispectrum and trispectrum can also be calculated, 
as described after Eq. (UH). 

Going to a super-large box, none of this may work. 
As we noted earlier, 5a might be big which would invali- 
date the calculation of its stochastic properties. Also, 5N 
might be big and then Eq. ([T7)) would apply only to the 
average of N within the box ((Tend then being a spatial 
average) , which might have little to do with the situation 
in the observable Universe. 



B. With o a curvaton-type fleld 

Now suppose instead that a has nothing to do with the 
inflation dynamics, as is typically (though not inevitably 
[T3 |) the case in the curvaton model. ^ Then at last we 
encounter a case where it may be useful to consider a 
super-large box, using what is often called the stochastic 
approach to the evolution of perturbations . 

In the stochastic approach one takes spacetime to be 
unperturbed with constant H (de Sitter spacetime). One 
smooths cr on a practically-fixed scale (1 -I- h)H with < 
& <C 1 a constant, and considers the probability Fit, a)da 
that cr(x, t) lies within a given interval. It satisfies the 
Fokker-Planck equation 



dF _ V dF 



d^F 



(23) 



^ For slow-roll inflation this corresponds to V' being much less 
than the corresponding quantity for the inflaton but we are not 
making any assumption about the inflation model. 
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This equation, applying to any slow-rolling field, corre- 
sponds to the Langevin equation describing the classical 
evolution plus the random walk zLH/2tt per Hubble time 
coming from the creation of the classical perturbation 
from the vacuum fluctuation. 

The point now is that the probability distribution may 
lose memory of the initial condition. In particular, if 
H/H^ is sufficiently small, F will settle down to 



F — const cxp — 



3H^ 



Via) 



(24) 



More generally, one can handle a significant variation of 
H within a given inflation model (see for instance [TtI . 

El). 

Reverting now to a minimal box, the probability dis- 
tribution F should apply to if the initial epoch is taken 
to be soon after the shortest cosmological scale leaves the 
horizon, since the box size is then not too many e-folds 
bigger than the Hubble scale at the initial epoch. In any 
case, one could calculate the probability distribution of 
a by going back to the Langevin equation. The very sim- 
plest case arises if tr is a pseudo Nambu-Goldstone boson 
(PNBG) with V negligible. Then it is deflned only in 
some interval < cr < /, and the noise term gives a an 
equal probability of being anywhere in the interval. 

Finally, given a probability distribution for a one may 
suppose that the actual value for a minimal box around 
the observable Universe is not too far from the most 
probable value (or from say the mean-square if we are 
not dealing with a PNGB). Of course this final step is 
speculative and may be modified by environmental con- 
siderations. Still, it seems that in this case the use of a 
super-large box purely to get a handle on the likely value 
of a within a minimal box may be helpful. 



V. CALCULATING THE CORRELATORS OF C 



Using the convolution theorem, 



(25) 



Eq. ^ determines the correlators of ( in terms of those 
of Sa. There is a sum of terms and the calculation is best 
done using Feynman-likc graphs [1, [l^, [l^l • A complete 
set of rules for constructing the graphs is given in [2Q| . A 
graph with n loops involves an integration of n momenta, 
while a tree-level involves no integration. 

Here I recall the estimates of the correlators made in 
d, Q for the quadratic truncation of Eq. ^ . To get an 
idea of what happens with higher terms included, I then 
consider the cubic truncation. In both cases I take da to 
be gaussian. (See for a loop contribution involving 
the bispectrum of Sa.) 



A. Quadratic truncation 

Truncating the field expansion after the quadratic term 
we have 



C(x) = Ar',5a(x) + iiV"5a2(x) 



(26) 



There are tree-level and one-loop contributions to the 
correlators of C,: 



(27) 
(28) 



^troc ^ 2N'^N"P„{ki)P^{k2)+ cyclic (29) 



loop 



(2^ 



/ d^pPaip)PaiPl)PAP2) (30) 



j,trcc ^ Ar'27v"2p^(fc^)p^(fc2)P,(fci4) + 23perms.(31) 



F 



loop 



1 N" 



8 (27r)3 
23perms. 



d''pP„{p)P, {p^)P, {p2)Pa {P2a) 

(32) 



We have defined pi = \p - ki|, p2 = \p + k2|, P24 = 
|p -I- k24| and ku = |ki -|- k4|. 

The 24 terms in Eq. ([31]) are actually 12 pairs of iden- 
tical terms, and the 24 terms in Eq. are actually 3 
octuplets of identical terms. The tree-level contribution 
to the bispectrum was given in [22| and the tree-level 
contribution to the trispectrum was given in Q (see also 
[2l|). The loop contributions to the spectrum, bispec- 
trum and trispectrum were given in Q , using Va instead 
of Prr a nd with Sa normalized to make N" = 1.^ (See 
also [2J| for the loop contribution to the spectrum of the 
axion isocurvature perturbation, given by an identical 
formula.) 

The subscript on the integral reminds us that Pa is set 
equal to zero at < L~^, cutting out a sphere around 
each of the singularities. This is necessary, because with 
Va perfectly flat there is a logarithmic divergence when- 
ever the argument of P^ goes to zero, ie. in the infrared. 
Allowing for slight scale dependence ofVa, infrared con- 
vergence is slow if it occurs at all. In contrast, there is 
good convergence in the ultra-violet for any reasonable 
behaviour of Va, and the integral will be insensitive to 
the actual cutoff fcmax- 

It is convenient to define what one might call a reduced 
bispectrum /nl and a reduced trispectrum by 



Be = ^hL[Pdki)Pc{k2) + cyclic] 



(33) 



Tc = -r^Fc(fci)^c(fc2)Pc(fci4) + 23perms.. (34) 



^ The SN formula was not invoked there, and indeed is irrelevant 
in the present context. All we really need is that f is some 
quadratic function of a gaussian quantity 5<p with zero mean. 
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At tree-level the reduced quantities are momentum- 
independent: 



j*tree 



6 
5 



" ftrcc 



(35) 
(36) 



In first order perturbation theory, this definition of /nl 
coincides with the original one [22] in first-order cosmo- 
logical perturbation theory, where /nl was defined with 
respect to the Bardeen potential which was taken to be 
<I> = |C. In that reference it was supposed to be inde- 
pendent of the momenta. Following |23| we will allow 
momentum dependence, and make the definition with- 
out invoking first-order cosmological perturbation the- 
ory."* The quantity here denoted as tq was introduced in 
and denoted as tnl. Taking them to be momentum- 
independent observation gives bounds |/nl| ^ 100 and 
tq < 10"*, which on cosmological scales are not expected 
to alter much if there is momentum dependence. 

Now I consider estimates of the loop contributions, tak- 
ing Va to be perfectly fiat and considering only the phys- 
ical regime ki 3> L^^ . One can arrive at an estimate of 
.pioop assuming that the integration is dominated by 
spheres around each of the two singularities, with radii 
of order k. This gives Q 



order as (7^^°°^ /"Pi"®). Special configurations of the mo- 
menta will give additional factors, but observation can 
probe only a fairly limited range of momenta. Pending 
further investigation of that issue, we conclude that the 
loop contributions are very suppressed for a minimal box, 
and unlikely to be observable. 

Our general finding that the loop contributions are 
negligible relies on the observational bounds on non- 
gaussianity, and holds because we assumed that Sa is 
the only field perturbation contributing to C,. As we have 
seen though, in the particular case that a is the infla- 
ton in a single-field slow-roll infiation model, the loop 
contribution is very small by virtue of slow-roll, and in 
consequence the non-gaussianity is very small. With this 
in mind, we can suppose that the infiaton field per- 
turbation gives the dominant contribution to C and write 



C = Cinf + N'5a{^) + -N"5a\^) + 



(39) 



where 5a is say the curvaton. Then the spectrum of Q is 
dominated by that of the first term, but its bispectrum 
and trispectrum might instead be dominated by the third 
term which would mean that they were dominated by a 
loop contribution In that case Ci„f could account for 
up to 90% of the total, without violating observational 
bounds on non-gaussianity. 



-jloop 



2N"^Vt\n{kL). 



(37) 



Cubic truncation 



We deduce that 



V. 



loop 



N" 



•^ticc 



VcHkL) <10-^HkL), (38) 



where the inequality comes from the observed spectrum 
T'f = (5 X 10^^)^ and the observational bound on /nl or 

(by coincidence those bounds give a similar result). 

The integration for 'p^°°P can actually be done ana- 
lytically [131 and it happens to give exactly this result. 
The loop integrals for the bispectrum and trispectrum 
cannot be done analytically, but they can be estimated 
in the same way as for 7?^°°p. Although a more general 
estimate could be made, we will take all ki to be of or- 
der a common value k for the bispectrum, and all ki and 
kij to be of order a common value for the trispectrum. 
Then, estimating the loop integrals from the contribu- 
tions of the singularities as we did for 'p^°°P^ on finds 
that (/n^Z/nD and {t°°^ /t^"''') are both of the same 



* We have no need of perturbation theory before horizon entry 
but it is needed to evolve perturbations afterward. Second order 
perturbation theory will be needed if |/nlI ^ li and S't that 
order <& and f are completely different functions. As a result 
/nl defined with respect to $ has nothing to do with the /nl of 
the present paper. Unfortunately, both definitions are in use. 



Now we include the cubic term: 

C(x) = N'Saix.) + -N"6cr^{x.) + -N"'Sa^x.). 
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(40) 



There are no additional tree-level diagrams for the spec- 
trum and bispectrum, but for the trispectrum [l^ the 
cubic truncation gives an additional term 



N"'N"'PAk2)Pa{k3)Pa{k4) + 3 perms. (41) 



The new term cannot cancel the old one since its de- 
pendence on the momenta is quite different. To extract 
optimal information from the observations one should re- 
define and introduce a new quantity (7^ by writing 



Tr 



1 



T^Pciki)Pcik2)Pciki4) + 23 perms. 



gc [PAk2)Pa{k3)Pa{k4) + 3 pcrms] . (42) 



(This is the same as the (7nl of |13| up to a numer- 
ical factor.) As with tq it is helpful to allow to be 
momentum-dependent because it can then correspond to 
a (slightly) momentum-dependent loop contribution. 
Pending the appropriate observational analysis, we will 



take ki 
Then 



k and keep the original definition of r^. 



^trcc3 



N'" 

77^ 



\n(kL) < W^, 



(43) 
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where the final inequahty assumes a minimal box and 
uses the observational bound on which we take to be 
the same as if were momentum-independent. 



Denoting the old contribution by ri"^""^ we have the 



ratios 



/•tree trcc2 trec3 
/NL ■ " C • " C 



N" 



N" 



N' 



N" 



iV'3- 



As pointed out in [T^l, t^'^'^"^ might be the first signal 
of non-gaussianity. This could happen in the curvaton 
model if the curvaton field evolves strongly after inflation. 

It could also in principle happen if a is the inflaton 
in a single-component slow-roll inflation model. Slow- 
roll requires only that all three slow-roll parameters are 
<^ 1, and one might have over a limited range of scales 
I^P S> |?7| and 3> e. In such a case though, we have to 
remember that the small non-gaussianity of 6a at hori- 
zon exit will be comparable with the non-gaussianity that 
we are considering here (ie. that generated by the non- 
linearity of the SN formula). The known estimates of 
the bispectrum [Tl| and trispectrum [l^ of da at horizon 
exit assume that is negligible, and will require mod- 
iflcation if it is not. Note also that in such a case the 
usual [25j formula n — 1 = 2r] — 6e the spectral index 
may fail as well (26j . Then one would have to rethink 
the the implication of the current measurement of n — 1 
and of the current bound on r = 16e, for non-gaussianity 
in slow-roll inflation. Of course, such a rethink is hardly 
going to alter the conclusion that non-gaussianity in this 
model will be very hard, if not impossible jl^l to detect. 

Now we turn to the loop contributions. In the presence 
of cubic and higher terms, one finds integrations over a 
single momentum. In the graphical representation, these 
correspond to loops which start and finish at the same 
vertex. It has been shown [20*1 that, instead of including 
such loops, one can replace the factors N', N" and so on 
by 



N' = N'{a) (iV'(x)) = N\ 



(45) 



and so on.^ When these 'renormalized vertices' are used, 
one need only draw 'renormalized graphs', which omit all 
lines corresponding to an integration over a single mo- 
mentum. We are working at cubic order, which means 
that only N' gets renormalized: 

iV'(x) = N' + N"6a{^)+'^N"'5a''{^) (46) 



N' 



N' + ^N"'{5a^), 



(47) 



with 



{6a^) = 



P,(fc)— ~7'.ln(fc,naxi). (48) 
L-1 k 



Let us verify that the loop contributions to the spec- 
trum are still suppressed at cubic order. The renormal- 
ized tree-level contribution is 



■tree tC^'S 



Using Eq. (j48|) and the bound Eq. ([43]) we find 



(49) 



(50) 



With the cubic truncation we have 1- and 2-loop con- 
tributions, which have no renormalization. The 1-loop 
contribution to V(; given by Eq. ((28|) . The 2-loop contri- 
bution is M 



2 — loop 



1 N'"'^ f 

6 (2^ A2P(gi)P(|q2-qi|)P(|k-q2|). 



.(51) 

Taking the integral to be dominated by the three singu- 
larities, we estimate 



2 — loop 



c 



j^_/pPa{p)PA\v-^) (52) 
]#]^^C^r'°°^ (53) 



N"'Vq 



VI. RUNNING 

We have advocated the use of a minimal box, but we 
did see that it might be useful to consider also a super- 
large box in order to get a handle on ct within a minimal 
box. Suppose that for some reason we decide to perform 
the whole calculation of the correlators in some super- 
large box with size L. We may compare the outcome 
of such a calculation with one done in some smaller box 
with size M <^ L, placed within the super-large box. (I 
am thinking of the size M as being minimal but that 
is not essential.) An interesting situation then arises, 
which was explored for the quadratic case in 0, [1].^ I 
now repeat that analysis in a different way, arriving at 
a differential equation instead of a finite-difference one, 
and extend it to the cubic truncation. Then I ask how 
the calculation may be of practical importance. 



A. General situation 

I think of the super-large box as having a fixed size 
L. For a calculation within the a smaller box of given 



5 The authors of 0| use A^' to denote N'{x), so that our A^' is 
equal to their (A'^'). 



In the labels L and M are interchanged so that L < M. In Q 
the labels are as here, but in Eq. (24) of Q it should be log{kM) 
instead of log(fcL). 
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size M and a given location, a of the previous equations 
becomes o'j\/, and Eqs. ^ and ([5]) become 



C(x) 



where 



and 



Sani'x.) = fT(x) - (7M, 



TV" 



(55) 



(56) 



(57) 



and so on. In terms of the original quantities we have 



SaM = {(7 - (tm) + Sa{x) 



(58) 



N'm = N' + N"iaM-W) + -N"'iaM~^f + --m 

For a calculation within the smaller box, correlators are 
defined (in position space) as averages over the smaller 
box. The spectrum of Sa is not affected, because SaM 
differs from 5a only by a constant. If Sa is gaussian, 
we can therefore forget about the change in box size as 
far as its stochastic properties are concerned. The same 
is not true of the correlators of ( though; they will be 
different because C(x) is a different function, and because 
the average is taken over a smaller region. Let us denote 
the spectrum defined within the smaller box by Pmc ^-nd 
similarly for the bispectrum and higher correlators. 

If we fixed the size and location of the smaller box 
within the original box, that would be the end of the 
story. But if the smaller box surrounds the observable 
Universe, it may be reasonable to suppose that we occupy 
a typical position within the original box. In that case, 
instead of considering Pm^ and so on, one might consider 
{Pmc) and so on, the quantities obtained by averaging 
over the location of the smaller box while keeping its size 
M fixed. One might hope that this average will give a 
reasonable estimate of the correlators, evaluated within 
a smaller box of size M that is fixed at our unknown 
location. 

Since the correlators calculated within a given box can 
be defined as spatial averages within that box, {Pmc) 
and so on must be equal to the quantities P^- and so on, 
evaluated directly within the super-large box. However, 
if Pmc and so on are calculated from Eq. ([55)1 and the 
spatial average within the super-large is then taken, the 
separation into tree-level and loop contributions is differ- 
ent. The loop contributions will increase with M, and 
the tree-level contributions will fall to compensate. As 
we noticed earlier, the tree-level contribution will usually 
dominate if the size M is minimal, but that need not re- 
main the case as M is increased to eventually become 
equal to the super-large box size L. In the following sec- 
tions we see how the compensation occurs, first for the 
quadratic truncation and then for the cubic truncation. 

The cosmological situation that we have described is 
analogous to one that occurs in quantum field theory. 



There, one also calculates correlators (usually time-order, 
corresponding to scattering amplitudes) which are the 
sum of a tree-level and loop contribution. To do the cal- 
culation one has to specify a renormalization scale Q. 
The correlators are independent of Q but the separation 
into tree-level and loop contributions is not. By choosing 
Q to be of the same order as the relevant energy scale (set 
say by the momenta in a scattering process) the tree-level 
contribution will normally dominate if it is present, oth- 
erwise the one-loop contribution will normally dominate 
and so on. The cosmological situation that we consider is 
similar, with Q replaced by M . The maximal value L of 
the super-large box, determined by the amount of slow- 
roll inflation long before the observable Universe leaves 
the horizon, provides the infrared cutoff of the theory. 
Its field theory analogue is the ultra-violet cutoff of the 
effective field theory, that dictates the maximum choice 
of the renormalization scale Q. 



B. Quadratic truncation 

Because Eq. ([55)1 is only quadratic, N" is just a number 
and the expectation value with respect to the super-large 
box is required only for the tree-level terms. We have 



c / 

ploop 

c 



(60) 



/ d^pP^{p)P„{\p-k\) (61) 

(62) 



(2^)- JM- 

(B*'™) = 2(iV;2)iV"P<,(fci)P4fc2)+ cyclic 



loop 



(2^)3 



d'pPa{p)Pa{pi)Pa{P2) 



(63) 



Be 

{Tl'--) = (iV;j)7V"^P^(A:i)Pc(fc2)n(fci4) + 23pern(fi4) 

1 iV"4 



12 \ AT "2^: 



T, 



loop 



c 



, d^pPaip)Pa{pi)Pa{p2)PAP2i) 
J M-^ 

(65) 



8 {2^Y jj^,-^ 
23perms 



Using Eq. we find 

{N'i,)^N'^ + N"^{{aM~^?). 



(66) 



Considered as a function of the position of the box with 
size M , aM — is simply (5ct(x) smoothed with a top-hat 
window function. Its mean-square is therefore 



{{-^M-^f) 



M-^ 



Paik) 



dk 

This expression generates M-dependence 
d 



d\nM 



{{aM-^?) = -r„[M-^). 



(67) 



(68) 



The M-dependence of the tree-level contribution to the 
spectrum is therefore 



J ptree 

dlnM 



iV"2p^(M-i)P,(fc), 



(69) 
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and similar expressions hold for the other tree-level con- 
tributions. 

Now we come to the loop contributions. Because phys- 
ical momenta have k ^ L^^, we can set P{p) — P{L^^) 
near a singularity of the integral ai p — 0. Differentiating 
the integral with respect to L, we see that the required 
cancellation occurs: 



d\nM ^ 



p 

d\nM < 



.loop 



(70) 



Similarly, taking into account all of the singularities, one 
can see that the same is true for the bispectrum and 
spectrum. 

Returning to the analogy with quantum field theory, 
Eq. ([55)1 is like the 'running' of a coupling constant (or 
other parameter) with the renormalization scale Q. This 
running makes the correlators of the field theory inde- 
pendent of Q. 



The running of the prefactor of the renormalized one- 
loop contribution Plj^°°^ is given by Eq. ([55)1 as 



dlnM 



(79) 



Taking into account the running of the integral calculated 
in the quadratic case this gives 



-f N"P,{M-^)P,{k). (80) 
Finally, the running of integral in p^-ioop gj^gg 

^^r^'w = \^,N'''MM-')Pi-^'^^-{k). 

(81) 

We see that the total running of vanishes as required. 



C. Cubic truncation 

I will just consider the spectrum. Evaluating it in the 
smaller box and then taking the expectation value with 
respect to the original box, Eqs. (l49l) - ((5T|) become 



P, 



P'r = {N'^i)PAk) 

1 (iV"2 



(71) 



1 — loop 



, d'pP.ip)PA\p-M) (72) 
2 (27r)^ J^,f-i 



P^ = - ,^ / d^qid^q2 



6 (27r)6 



X P(-7i)P(|q2-qi|)P(|k-q2|). (73) 

The renormalized vertex in the smaller box is 

1 



{Sa^),, - / P^{k)^ 



dk 
k 



(74) 
(75) 



Using Eq. ([59]) this gives 



(iV;2) = N'^+{N"^ + N'N"') {(aM-^f) + \N"'^(Sa^f 

(76) 

Only the middle term is M-dependent, giving the run- 
ning 



(iV;|) = 7V"2p,(Af-i). 



(77) 



dlnM 

The running of the renormalized tree-level contribution 
Pj5j^° is therefore simply 

-^PmT = ~N"P4M-')P4k), (78) 
the same as in the quadratic case. 



D. Application 

Does the running have a useful application? At first 
sight the answer would seem to be 'yes', because by go- 
ing down to a minimal box the loop contributions become 
negligible. Unfortunately, the gain is illusory because we 
are not actually calculating correlators within any par- 
ticular minimal box. Instead we are calculation the ex- 
pectation values of the correlators within a minimal box 
(taken within the super-large box). But these are just 
the actual correlators calculated within the super-large 
box. As a result, the calculation has all of the uncertain- 
ties, and possibly fatal problems, that come with the use 
of a super-large box. 

The problem is that the correlators calculated within 
a super-large box may be quite different from the ones 
observed in our Universe. If we throw down a minimal 
box within the super-large one, the correlators calcu- 
lated within the minimal box will depend on its location. 
There is no reason to think that a particular correlator, 
calculated with the minimal box at our location, will be 
very close to the result obtained by averaging the posi- 
tion of the minimal box. There is even less reason to 
think that such will be the case simultaneously for all 
correlators. 

To quantify this concern one would like an estimate of 
the likely difference between the averaged correlator and 
the one observed. Extending the terminology coined for 
the cmb multipoles, one may call that cosmic variance. 
It will be defined by the correlators evaluated within the 
super-large box. 

As an example we may consider the simplest curva- 
ton model, where the curvaton has a quadratic potential 
and H/H^ is negligible (see [5] for this case and further 
references).^ In that case, W vanishes in a super-large 



Taken literally this case is not realistic because it gives spec- 
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box and um is generated entirely by the perturbation 
5a Ml having the gaussian probabihty distribution given 
by Eq. ([24| . According to taste, one may fold in this a 
priori expectation with environmental considerations. 



VII. CONCLUSION 

I have explained how the use of a minimal box leads 
to fairly clean predictions. In particular I have verified 
that it makes some specific loop contributions small, by 
virtue of observational constraints on non-gaussianity. 

I have also pointed to some of the uncertainties and 

tral index bigger than 1. Keeping H / negligible, the required 
value n = 0.95 could be generated by giving the quadratic po- 
tential a bump in the middle. Alternatively one could invoke sig- 



possibly fatal problems, that may come with the use of a 
super-large box. Some of the issues raised here are quite 
deep and more work needs to be done. Provisionally 
though, it would seem that the only use of a super-large 
box is in its possible provision of a probability distribu- 
tion, for the average of a curvaton-type field within a 
minimal box. 
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nificant H / though the probability distribution of ctjv/ would 
not then be given by Eq. II 241 1. 
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